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MODELS FOR THE ANALYSIS OF CN VIOLET RADIATION BEHIND 
SHOCK WAVES IN AIR CONTAMINATED WITH 
CARBON-BEARING COMPOUNDS 

By John E. Nealy 
Langley Research Center 

SUMMARY 

Techniques for analyzing the effects of contamination radiation from the CN violet 
band system behind normal shock waves are presented. The problem is approached in 
two ways: (1) a simple chemical system is set up for a known contaminant, and (2) a 

computer program which utilizes a free-energy minimization technique is used to eval- 
uate the CN particle concentration for given atom percentages of carbon. Equilibrium 
radiation from CN and air with various contamination levels is determined for shock 
velocities ranging from 20 000 to 44 000 ft/sec (6.10 to 13.41 mm/jusec) at free-stream 
pressures of 0. 1 and 1 mm Hg. Radiation from CN in the chemical relaxation zone behind 
a normal shock with C 2 N 2 as the free-stream contaminant is also considered. Results 
of the calculations are presented graphically and discussed. 

The results indicate that radiation from contaminated air in the velocity range 
between 20 000 and 30 000 ft/sec (6.10 and 9.14 mm/qsec) may be dominated by small 
contamination levels. However, at speeds greater than about 33 000 ft/sec 
10.06 mm/jnsec), relatively high contamination levels can be tolerated for equilibrium 
lows because the contamination radiation is minor in comparison with air radiation. The 
alculations for the nonequilibrium condition further indicate that the ratio of CN radia- 
ion to radiation from air in the chemical relaxation zone is comparable to the ratio for 
he equilibrium condition, and that contamination radiation probably presents no greater 
iroblem in the nonequilibrium region than in the corresponding equilibrium region. 

INTRODUCTION 

Carbon -bearing compounds in the presence of nitrogen at elevated temperatures are 
generally involved in a rather complex chemistry. For temperatures greater than about 
3000° K, one of the byproducts is the cyano radical CN, which gives rise to a relatively 
intense radiation in the near-ultraviolet and visible -violet region of the spectrum. Very 
often, conditions are such that the carbon-bearing contaminants are difficult to remove 



from the system, and the CN violet radiation becomes a severe hindrance to spectro- 
graphic studies (refs. 1,2, and 3). 

The specific purposes of this report are to investigate the behavior of the CN con- 
tamination radiation for a fairly wide range of working conditions, and to provide an aid 
in defining tolerable free-stream contamination levels for radiation heat -transfer studies 
in shock tubes and expansion tubes. The problem of fulfilling these purposes has been 
approached in two ways. In the first method, the production of CN behind a normal shock 
is assumed to arise from small amounts of cyanogen C 2 N 2 in the free stream. Calcula- 
tions for both equilibrium and nonequilibrium conditions have been made for this assump- 
tion. A second method allows the shocked gas to contain a specified ratio of carbon atoms 
to air atoms and uses the minimum of free energy to define the resulting chemical equi- 
librium. The results of these two types of calculations are presented for various 
amounts of contaminant and over a range of shock speeds from 20 000 to 44 000 ft/sec 
(6.10 to 13.41 mm/psec). These results should be an aid in analyzing observed contami- 
nation radiation from CN and will be indicative of the effects to be expected from a given 
contamination level. 


SYMBOLS 


F 

h 

j 

kf> k b 
K 1 

K 2 

L 

M 

n 


radiative flux, W/cm 2 

Planck constant, 6.63 x 10“^ J-sec 

emission coefficient, W/cm 3 -4ir sr 

forward and reverse chemical rate constants, cm 3 /mole-sec 

n 2 

11^1 -vj 

chemical equilibrium constant, , moles/cm 2 

n C 2 N 2 

Upl Ili^r _ 

chemical equilibrium constant, — , moles/cm 3 


CN 


shock layer thickness, cm 


constant, defined by equation (1), moles/g 
molar concentration, moles/cm 3 


pressure, mm Hg 


2 



r nose radius, cm 

t time, sec 

T temperature, °K 

V g shock speed, ft/sec (mm/psec) 

x mole fraction 

6 shock detachment distance, cm 

PC number of carbon atoms per molecule 

number of unspecified atoms per molecule 
v frequency of photon, Hz 

£,y distance behind shock front, cm 

p density, g/cm^ 

4> atom fraction of carbon 

Subscripts: 

1 free-stream conditions 

2 conditions behind shock wave 

o standard atmospheric conditions 

eq equilibrium conditions 

An asterisk (*) indicates an electronically excited state. 
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CALCULATIONS AND RESULTS 


Determination of Radiance of CN Violet Band System 

For a given number of CN particles whose excited states are populated according 
to a known gas temperature, the mean radiation per molecule for a particular transition 
may be determined from the oscillator strength, or f-number, of the transition along with 
the spectroscopic constants pertaining to the levels involved. For temperatures above 
approximately 4500° K, the CN violet band system resulting from the transition 
B 2 2 - X 2 2 is the predominant radiation from this molecule, with major band progres- 
sions located at 4216 A, 3883 A, and 3590 A. The calculation of the radiance per mole- 
cule as a function of temperature for this band system was performed by Henry T. 
Woodward at Ames Research Center, using a computer program similar to that described 
in reference 4. This type of calculation is commonly referred to as the "smeared rota- 
tional line" model, and has been shown to be a very good approximation to the more 
refined single-line calculation when the band system is reasonably compact (ref. 5). The 
radiance of CN violet used in this paper is graphically presented in figure 1. It corre- 
sponds to an f-number of 0.020 and a heat of formation of CN of 8.2 eV. 

Chemistry of CN Behind Normal Shock Waves 

In certain radiation experiments performed in the Ames combustion-driven expan- 
sion tube, high levels of CN violet radiation were observed (ref. 1). This radiation 
occurred when mylar and teflon diaphragms were employed and indicated that foreign 
substances from the diaphragms containing carbon appeared in the test gas. In one par- 
ticular case when air was used as the test gas, the total gas -cap radiation was measured 
as 23 watts for a model in a flow of speed 7.33 mm/psec at Pj/p Q = 4.1 x 10“ 3 . For a 
ballistic range test under similar conditions (V s = 7.52 mm/p sec at P\/ P Q = 4.1 x 10- 3 ), 
gas-cap radiation was observed to be 2.4 watts. This value is essentially a decade below 
that obtained in the expansion tube. From spectral analysis of the expansion-tube radi- 
ation, one -third of the integrated spectrum was estimated to be due to the CN violet band 
system, while most of the remainder consisted of atomic line radiation from metallic ele- 
ments. Under these test conditions, the mole fraction of CN required to produce the radi- 
ation attributed to the violet band system is 1.8 x 10" 3 . 

Since a certain amount of contamination in the free stream may be unavoidable, it 
is of interest to learn the effects of CN for a range of operating conditions. One approach 
to the analysis is to allow the formation of CN behind a normal shock to result from a 
known amount of C 2 N 2 contamination in the free stream. It is further assumed that the 
state variables of the shocked gas attain the same equilibrium values as for pure air. 
Normal-shock tables were used to determine these values (ref. 6). For a first 
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approximation to the chemistry of the contaminant, the reactions considered were merely 
the formation and dissociation of CN: 


K 1 

C 2 N 2 ^ 2CN 
k 2 

CN ^ C + N 


The equilibrium constants, defined as 


Kl 


n 


2 

n CN 

c 2 n 2 


and 


k 2 


n C n N 

n CN 


are functions of temperature only and may be computed from partition functions or tabu- 
lations of the free energy. The equilibrium constants, along with the conservation rela- 
tion for carbon atoms 


\ ( 2n C 2 N 2 + n CN + n d = Constant s M 


( 1 ) 


may be used to evaluate the CN concentration at a given temperature and density as 
follows: 


n CN" 


-( k 2 + n N ) + ^( K 2 + u n ) 


8pMn^ 

+ ~*T 



, pMn N (m«j\ 

K 2 + n N \Kl / 

The results of this approach are illustrated in figures 2 and 3 for free-stream mole 
fractions of C 2 N 2 of 0.0001, 0.001, and 0.005 for values of Pj of 0.1 and 1 mm Hg with 
a free-stream temperature of 300° K. The concentration of nitrogen atoms behind the 
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shock was assumed to have its pure-air value and was determined from equilibrium 
chemical properties of high-temperature air (ref. 7). Values of air radiation used in the 
figures were taken from Allen (ref. 8). Allen does not give the emission coefficient 
directly but plots the total radiative flux from one side of an infinite air slab of 1-cm 
thickness, taking absorption effects into account. For an optically thin gas, the flux F 
from an infinite slab is given by 



where L is the slab thickness. Since this relationship was used in obtaining the emis- 
sion coefficient for air in figures 2 and 3, some inaccuracy exists in the air curves at 
conditions for which absorption effects are significant. These assumptions and approxi- 
mations overestimate the CN concentration because the reactions involving CO, O, and the 
ionization of C have been neglected. Such an approach is admittedly only approximate, 
but this method of calculation is extremely simple and does show the trend to be expected 
at higher velocities, namely that the CN radiation becomes negligible as V s increases. 

The problem was approached in a different way by utilizing a computer program 
which determines the chemistry behind the shock by minimizing the free energy of the 
system. It was assumed that a small percentage of an unknown carbon-bearing compound 
was present in the free stream but had no effect on the cold-gas properties. A value was 
given to the ratio of carbon atoms to the total number of atoms in the gas, and the reaction 
chemistry of the carbon behind the shock was determined according to the free-energy 
minimization process. This calculation procedure is physically equivalent to a model in 
which the carbon-bearing contaminant is dissociated into its atoms with no energy drain 
from the shock wave and in which only the carbon atoms are chemically active with the 
ordinary air constituents. Again, this model would tend to provide an upper limit to the 
CN concentration, since intermediate species involved in breaking up the contaminant have 
not been considered. An advantage to this approach is that some idea of an acceptable 
contamination level for a carbon-bearing compound may be obtained' from the results of 
the calculations. For example, if a contaminant A contains Hq carbon atoms per 
molecule and jU u atoms of additional unspecified elements, the atom fraction of carbon 
is 


M C X A 

(M u + Mc) x A + 2 ( x N 2 + x 0 2 ) 
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with 


X A + x N 2 + X C>2 = 1 
If 

X A <<: x N 2 + x 0 2 


the mole fraction of the compound is approximately 



( 2 ) 


Calculations have been made for atom fractions of carbon of 0.001, 0.005, 0.01, and 0.05 
for pj values of 0.1 and 1 mm Hg. The results are shown graphically in figures 4 
and 5. In these figures, the CN radiation reaches a maximum for V s ~ 26 000 ft/sec 
(7.92 mm/psec) and decreases rapidly thereafter. This trend, along with the rapid 
monotonic increase of the air radiation, illustrates the fact that the contamination radi- 
ation becomes negligible at the higher speeds. 

For V s = 23 000 ft/sec (7.01 mm/p sec), the CN radiation observed in the Ames 
expansion tube is approximately that calculated for 0= 0.05, for equilibrium conditions 
behind the shock. (See fig. 5.) Equation (2) applied to a compound containing 10 carbon 
atoms per molecule gives a contaminant mole fraction of 0.01. Thus only 1 -percent con- 
tamination of the test gas by such a compound could conceivably produce behind a shock 
CN violet radiation of an intensity several times that of air. However, radiation caused 
by this species at speeds greater than about 32 000 ft/sec (9.75 mm/psec) could be 
completely ignored as far as radiation heat transfer is concerned. At the higher speeds, 
atomic and ion continuum radiation from the dissociated constituents of a contaminant is 
not expected to be significantly different from radiation from atoms and ions of pure air, 
since theoretical relationships for the various types of continua predict a weak depend- 
ence of the emission coefficient on the atomic weight of the species. 

The CN violet radiation as calculated by the previously described methods is com- 
pared to a direct determination of in which the equilibrium mole fraction of CN was 

computed by using the normal-shock, free-energy minimization program of reference 9 
(figs. 6 and 7). This program was applied to air with a free-stream mole fraction of 
C 2 N 2 of 0.005. These results are compared with the first approximate solution for 
( X C 2 N 2) = 0.005 and with the solution for a corresponding atom fraction of carbon of 

0.005 at free-stream pressures of 0.1 and 1 mm Hg. Both methods predict accurately 
the velocity regime in which the CN radiation begins to decrease. Below 26 000 ft/sec 
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(7.92 mm/fxsec), the first approximate solution leads to an overestimate of the number 
density of CN behind the shock, primarily because of neglect of the exchange reaction 
CN + 0 : CO + N. However, the results for <£= 0.005 agree very well with the direct 
calculation over the entire velocity regime for which CN radiation is most significant. 

Nonequilibrium Calculations 

The problem of nonequilibrium radiation from CN was considered also. A computer 
program which calculates the finite-rate chemistry behind a normal shock was used to 
compute the concentration of CN as a function of distance. The pressure and temperature 
were assumed to increase discontinuously across the shock as determined from the con- 
servation equations for mass, momentum, and energy. The reaction mechanism included 
the reactions considered most important for air, in addition to the following reactions 
involving the carbon compounds: 


C 2 N 2 + M £ 2CN + M 
CN + OsCO + N 
CO + O + Ma C0 2 + M 
CO+M=C+0+M 
CN + MsC + N + M 
C + e a C + + 2e 

Rates for reactions occurring in air were obtained from reference 10; the dissociation 
reaction used for C 2 N 2 was that given in reference 11. The remaining reaction rates 
were taken from reference 12, with the exception of the ionization of carbon, for which a 
classical electron-impact formula was used to define the ionization rate. 

For a first approximation to the radiation from the non equilibrium zone, the internal 
energy modes were assumed to be in equilibrium with the local translational temperature. 
In order to obtain a more realistic model for the nonequilibrium radiation, it is necessary 
to account for the finite excitation time of the excited electronic levels of the radiating 
species. In the case of radiation from air, the experimental data of reference 13 were 
used to determine the distance behind the shock at which the radiation intensity was max- 
imum. The intensity was then assumed to increase linearly with distance up to this point, 
and thereafter to be determined by the local translational temperature. 

2 

For the case of the CN violet, reaction rates involving the excitation of the B E 
level of CN and the deexcitation by collisions and by spontaneous emission were employed 
according to the method of reference 14. The profiles for the CN violet were calculated 
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by this method until the resulting intensity approached that obtained when the B^S level 
was populated according to the local translational temperature. The additional reactions 
and rates used for these calculations are 


CN + M s: CN* + M 


CN* — CN + hv 


7kf = (8.62 x io-I^T 5 - 5 ^ 100/T... cm 3 \ 
^ mole -sec/ 

Vk b = (8.62 x 10-1^5.5° — — ) 

^ u mole-sec/ 

(k f = 8.86 x 10 6 sec-l) 


For each set of initial conditions, the emission coefficients of air and CN violet and 
the temperature and mole fraction of CN, along with the equilibrium asymptotes, are 
illustrated as functions of the distance behind the shock (figs. 8 and 9). The equilibrium 
values for T and were determined from the equilibrium normal-shock program 

(ref. 9), and the equilibrium air-radiation values are those obtained from Allen (ref. 8). 
The nonequilibrium air radiation was calculated from correlation formulas which give the 
emission coefficient per particle as a function of temperature for various transitions. 
These formulas for the air species are given in reference 15, and the CN violet radiation 
was calculated from a correlation formula determined from the data in figure 1. Need- 
less to say, these formulas are extrapolated in certain instances to extremely high tem- 
peratures where they cannot be expected to remain very accurate. Furthermore, the 
oxygen and nitrogen line radiation and free-bound radiation, which are major contributors 
at high temperatures, have been omitted from the air-radiation calculations because 
acceptable correlation formulas were not immediately available. The omission of this 
radiation contributes to the variance between the value of the emission coefficient of air 
at 5 cm and the corresponding equilibrium value given by Allen for a velocity of 
33 000 ft/sec (10.06 mm//isec). (See fig. 8.) 

One conclusion which may be drawn from these calculations is that the CN violet 
contamination radiation presents about the same problem in the chemical relaxation zone 
as in the equilibrium zone, at least for the lower shock speeds. However, the CN number 
density profile may differ considerably for various contaminants, and it is expected that 
the CN violet profile will not exhibit as rapid a rise in the nonequilibrium zone when a 
contaminant does not contain the CN radical in its molecular structure. To illustrate 
this, the nonequilibrium calculation has been performed for 1 percent CO 2 contami- 
nation. The resulting profile for the CN violet radiation is shown in figure 10 for 
V s = 23 000 ft/sec (7.01 mm/psec) and p x = 1 mm Hg, and is compared with the calcu- 
lation for 0.5 percent C 2 N 2 contamination for the same speed and initial pressure. The 
equilibrium value shown is for (p= 0.005. 
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The nonequilibrium normal-shock calculations may be used to estimate the effects 
to be expected for a blunt body in expansion-tube flow by assuming a linear velocity pro- 
file through a gas cap with shock detachment distance 6. The particle residence time in 
the shock layer may then be equated to the particle time after passage through a one- 
dimensional normal-shock wave. Use of the empirical expression for shock detachment 
distance for a spherical -nose radius r 


P i 

6 ~ 0. 8 — r 

p 2 

together with the assumption of a linear velocity profile between the shock and the stag- 
nation point, leads to the following approximate expression for particle residence time in 
the shock layer at a distance £ behind the shock front: 

‘“-vH 1 -!) 

For a typical case of r = 3 cm, V s = 23 000 ft/sec (7.01 mm/psec), and £ = 0.756, the 
calculated particle time is 3.8 /isec. The distance behind the shock in the standing 
normal-shock calculation corresponding to this particle time is about 0.35 cm. Fig- 
ures 10, 11, and 12 show that for these conditions, nonequilibrium effects must be taken 
into consideration. In order to illustrate this further, the radiation profiles have been 
integrated, and the ratios of contaminant radiance to air radiance are compared in table I. 
The ratio for the corresponding equilibrium case is also given. 

TABLE I.- RATIOS OF CONTAMINANT RADIANCE 
TO AIR RADIANCE 


Contaminant 

p0.35 

J 0 Jcn d y 

f0.35 

Jq ^air dy 

0.005 C 2 N 2 

0.212 

.01 co 2 

.00452 

<p = .005 (equilibrium only) 

.187 
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The C 2 N 2 contamination level of 0.005 is seen to lead to a nonequilibrium flux ratio about 
the same as the corresponding equilibrium ratio, whereas the CO 2 contaminant leads to a 
much lower flux ratio. 


CONCLUDING REMARKS 

Techniques for analyzing the effects of contamination radiation from the CN violet 
band system behind normal shock waves are presented. Equilibrium radiation from CN 
and air with various contamination levels is determined for shock velocities ranging from 
20 000 to 44 000 ft/sec (6.10 to 13.41 mm/jnsec) at free-stream pressures of 0.1 and 
1 mm Hg. 

For all equilibrium cases, the CN radiation is seen to be negligible at the higher 
velocities, even for the highest contamination levels considered, because of the depletion 
of the CN particle concentration by dissociation. However, at the lower speeds, relatively 
small amounts of contaminant may cause the air radiation to be completely dominated by 
the CN violet band system. 

The primary conclusion drawn from the nonequilibrium calculations is that the CN 
violet contamination radiation can be about as severe in the chemical relaxation zone as 
for the corresponding equilibrium case. Of course, C 2 N 2 is an unrealistic shock-tube 
contaminant and was used in the calculations only as a convenient method of obtaining CN 
behind the shock. More realistic contaminants might include various forms of hydrocar- 
bon oils and perhaps graphitic carbon. For these types of contaminants, the production 
of CN behind the shock would not take place nearly as rapidly as for the calculated cases 
and, in many instances, the CN number density would be monotonically increasing to its 
equilibrium value. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., November 17, 1966, 

129-01-08-35-23. 
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CN violet emission coefficient, j- N , W/particle-4n sr 



Figure 1.- Average CN violet emission coefficient per particle for an f-number of 0.020. 
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Emission coefficient, j, W/cm 3 -4 



Figure 4.- 


















temperature and mole fraction of CN for initial contamination level of 0.5 percent C 2 N- 
/ s — 33 000 ft/sec U0.06 mm/Msec); p, = 1 mm Hg. 










Distance behind shock, cm 




